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qh' Abstract 

Males are often the "sicker" sex with male biased parasitism found in a 
►^ , taxonomically diverse range of species. There is considerable interest in the 

O ! processes that could underlie the evolution of sex-biased parasitism. Mating 

^ ' system differences along with differences in lifespan may play a key role. We 

examine whether these factors are likely to lead to male-biased parasitism 
[~>^ . through natural selection taking into account the critical role that ecological 

O \ feedbacks play in the evolution of defence. We use a host-parasite model 

with two-sexes and the techniques of adaptive dynamics to investigate how 
mating system and sexual differences in competitive ability and longevity 
can select for a bias in the rates of parasitism. Male-biased parasitism is 
selected for when males have a shorter average lifespan or when males are 
5^ \ subject to greater competition for resources. Male-biased parasitism evolves 

as a consequence of sexual differences in life history that produce a greater 
proportion of susceptible females than males and therefore reduce the cost of 
avoiding parasitism in males. Different mating systems such as monogamy, 
polygamy or polyandry did not produce a bias in parasitism through these 
ecological feedbacks but may accentuate an existing bias. 

Keywords: Life-history evolution, male-biased parasitism, adaptive 
dynamics, evolution of disease resistance 
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1. Introduction 



Across a diverse range of taxa tliere is often a sex bias i n parasitism rates, 



with males commonly the "sicker" sex (jZukl . Il990l . |2009| ). Although male 



biased parasitism is by far the mos t commonly reported s ex-biased parasitism 



(iPouliru . Il996l : ISchalk and Forbesl . 119971 : iPerkins et al.l . l2003l : iFerrari et al. 



20041 ). higher rates of parasiti sm have also been reported for females in avian 



hosts (JMcCurdy et al.l . Il998l ). It remains unclear what processes can ac- 



count for sex-biased parasitism and in particular the higher prevalence of 
disease in males. One possibility is that males exhibit different behaviour 
that leads to greater exposure (i.e. larger home ranges or more risk of 
i nfection for any given exposu r e bec ause of damage caused by fighting) 
( Bundyl . Il988l : iResif and Amosl . 120101 ). Bias may also result from under- 
l ying differences in life- h istory characteristics between males and females 
( JMoore and Wilsoru . |2002| ) including the idea that the larger physical size and 
growth ra tes of males make them a more accessible and attractive target for 
parasites (JMoore and Wilsonl . 120021 ) . There is also clear evidence of a physio- 
logical basis for differences in susceptibility with for example androgenic hor- 
mones i n males (testosterorie in vertebrates), acting to depress the immune 



syste m ( iMoore and Wilsonl . l2002l : iFolstad and Karterl . Il992 
1988h . 



Alexander and Stimsod . 



Beyond physiological mechanisms, it has been proposed that life-history 
theory could explain immune differences from an adaptive point of view in re- 
lation to sex-specific reproductive strategies. In particular it has been argued 
that the reduced investment in susceptibility is due to trade-offs between male 
mating effort and immune defense. In this scenario as the strength of sexual 
selection on males increases, the magnitude of the sex differences in immuno- 
competence will increase. In essence the argument is that a reduced immune 
syste n as may be t he unavoidable price of being male due to sexual selection 



( iZukl . Il990l . I2OO9I ). Following on from this a polygamous mating systems 



should lead to greater differences in male biased parasitism and that pe rhaps 
under polyandry females should be more susceptible (IZukl . Il990l |2009[ ). The 
ba sic assumptions of t hese i deas have recently been examined theoretically 
by IStoehr and Kokkd (120061 ) who determine optimal allocation of resources 
between immunity, survival and reproduction in males and females, under 
varying levels of sexual selection. This work has given a strong theoretical 
underp inning to ideas that sexual selectio n can explain male-biased para- 
sitism (IStoehr and Kokkd . l2006l : IZukl . |2009| ). In addition, iMoore and Wilson 



(120021 ) carried out a meta-analysis using two measures of the strength of 
sexual selection-mating system and sexual size dimorphism-that showed that 
sexual selection was associated with sex differences in parasitism. 

There is a large body of theoretical work that has emphasised the im- 
portance of ecological fe e dback s to the evolution of host defence to infectious 



disease (see lBoots et al.l ( 120091 ) for a review). It is clear from this theory that 
host life-history is critical to level of defence that evolves. In particular, it 
is often, although not always the case, that increased resistance to parasites 
is more likely to evolve for long-lived hosts ( iMiller et al.l . l2007al ). Differences 
between males and females in terms of their life-histories may therefore be 
enough to explain the evolution of different levels of investment in defence 
through natural rather than sexual selection. In particular differences in 
ecological feedbacks between males and females due to differences in their 
life-histories may underpin the evolution in reduced investment by males in 
defence and therefore result in higher transmission of infection. Furthermore, 
the ecological feedbacks due to monogamous and polygamous mating systems 
have different effects on the evolution of male and female investment in de- 
fence. In particular, different mating systems cause different patterns in the 
way in which densities of males and females feedback into the evolutionary 
dynamics. 



A recent model (JResif and Amosi . l2010l ) has shown the importance of epi- 
demiological feedbacks in determini ng male-biased parasitism through natu- 
ral selection. iResif and AmosI (120 lOl ) focus on trade-offs in the disease charac- 
teristics such that differences in exposure may lead to differential investment 
in resistance. We develop a continuous time host-parasite system that rep- 
resents two-sexes and investigates how the mating system and differences in 
competitive ability and longevity between sexes influences the level of re- 
sistance to infection that evolves. As such we are examining the effects of 
mating system and host life-history under natural rather than sexual selec- 
tion. Crucially we include epidemiological dynamics and focus therefore on 
how ecological feedbacks affect the evolutionary process. We find that differ- 
ences in lifespan are enough to explain male biased parasitism. Differences 
in the mating system act only to accentuate an existing bias. Our work fur- 
ther emphasises the importance of including epidemiological feedbacks when 
studying the evolution of defence. 



2. Methods 



The underlying host-parasite framework is based on the classical ap- 
proaches for modelling the population dynamics of directly transmitted mi- 
croparasites (see Andersonl981) and which hav e been successfully extended 
to understand the evolution of host resistance (JBoots and Haraguchil . Il999 : 



Boots et al.l. 120091) . This fra mew ork is modified follo wing the techniques of 



Lindstrom and Kokkol (119981 ) and lMiller et al.l (l2007bl ). to represent a two-sex 
host parasite model that considers the dynamics of males and females sepa- 
rately. This is achieved by representing births as the harmonic mean function 
proposed by lCaswell and Weeks! (119861 ). that depends on the densities of the 
two sexes and declines to zero in the absence of either sex. This function 
can also be modified to approximate different mating systems (monogamy, 
polygyny and polyandry). The theoretical framework is therefore represented 
by the following system of nonlinear ordinary differential equations for the 
densities of susceptible, S, and infected, /, males and females, represented 
by the subscripts m and / respectively. 



dbm 

IT 

dim 

~df 
dSj_ 

dt 

dif 

dt 



-B{Sm, Sf){l — QmH) — bmSm — PmSmi^f + ^n 

Pm>~'m\l f ~r Im) \'^m ~r Oijlm 

^BiSm, Sf)il - QfH) - hfSf - PfSfilf + I J 

^ = l3fSf{If + Im)-{bf + a)If 



(1) 



Where H = Sf + Sm + If + Im is the total host density. Births are 
divided equally between males and females according to the harmonic birth 
function, B{Sm, Sf) , which describes the dependency of the birth rate on the 
density of either sex and mating strategy. The birth rate is modified due to 
density- dependent competition for resources with the parameter q, and the 
population has a natural death rate, b. Infection can occur through contact 
between susceptible and infected individuals with transmission coefficient, 
/3, and the disease induces an additional mortality while infected at rate a. 
(The subscripts on some parameters distinguishes between male and female 
specific parameters.) The harmonic birth function, B{Sm,Sf) , is derived 



from (jCaswell and Weeksl . Il986l ). and takes the following form. 
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Here, Cm and Cf represent the contribution that males and females make 
to the birth rate and h represents harem size and can be manipulated to 
represent different mating systems. When h > 1 it represents a polygenic 
mating system (births are maximised when females exceed males), when 
h < 1 it represents polyandry (births are maximised when males exceed 
females) and when h = 1 it represents monog amy (births are maximise d 
when males and females are equally abundant) (jCaswell and Weeksl . Il986l ). 

To examine the ev olution of parasite resistance we follow the techniques 



of adaptive dynamics (JGeritz et al.l . Il998l : iBoots et al.l . 120091 ) . We assume a 



'mutant' strain of host can occur at low density and attempt to invade the 
established 'resident' strain which is at its equilibrium density. The mutant 
male host strain differs from the resident strain in terms of its transmission 
coefficient (3^ compared to Pm for the resident (a similar difference can occur 
for the female transmission coefficient and we will use " ~ " to represent 
the mutant parameters). In line with previous studies into the evolution 
of host resistance it is assumed and that a benefit in terms of increased 
resistance to infection is bou ght at a cost in terms of a reduced birth rate 
( JBoots and Haraguchil . Il999l ). For this study we impose the trade-off Cm = 
g{(3m) and c/ = g{(3f) ■ The trade-off is defined by 
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Which is a smooth curve between the minimum and maximum values of 
the birth and transmission rates and in which the parameter 7 controls the 
curvature (and therefore cost structure) of the trade-off. (We will discuss 
other possible trade-offs later.) 

The fitness is the long-term exponential growth rate of a phenotype 
in a given environment. We initially consider the situation where the fe- 
male parameters are fixed and we allow the male parameter (3^ (and Cm 
via the trade-off) to evolve. A proxy for the fitness, i?, can be calcu- 
lated as the determinant of the jacobian m atrix, J, at the st eady state 
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and therefore, R, can be represented by the following expression 
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The fitness proxy, R, can be used to determine the position of evolutionary 
singular points and the evolutionary behaviour at the singular point. Evolu- 
tionary singular points are determined when the fitness gradient ^^ = 

OPm a — « 

Hm — Mm 

which equates to solving the following expression 
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The evolutionary behaviour at the singular point is determined by analysing 
the second order p artial derivat i ves of R with respect to the mutant and res- 



ident parameters ( jGeritz et all Il998l ). Previous studies have assessed how 



the trade-off cost structure or underlying population dynamics can influence 
the evolutionary behaviour and also induce evolutionary branching leading t o 
diversity in host strategies ( iBoots and Haraguchil . Il999l : iBoots et al.l . 120091 ) . 
The focus here is to examine whether different levels of resistance can evolve 
between males and females and therefore to allow us to concentrate on this 
issue we ensure that the underlying population dynamics are point equilib- 
rium and that the trade-off has sufficiently accelerating costs that the singular 
point is an evolutionary stable attractor. We can then assess how the position 
of the singular point changes as other life history parameters are varied. 

To examine the coevolution of male and female resistance properties we 
determine the female fitness function (which depends on the evolving pa- 
rameter Pf and Cf via the trade-off). The male singular points are plotted 
against /3f (for a fixed female strategy using the method outlined above) and 
the female singular points (against a fixed male strategy) are plotted against 
I3m ■ The int ersection of these l ines p roduces a coevolutionary attracting sin- 
gular point ( JRestif and Koellal (120031 ) , note that again the trade-off is chosen 
to ensure the population dynamics exhibit a point equilibrium and that the 



singular point is a coevolutionary stable attractor). Using this method it 
is possible to determine how the coevolutionary singular point varies with 
changes in underlying life history parameters. 

3. Results 

The evolutionary behaviour is dependent on feedbacks that arise in the 
ecological dynamics and therefore it is useful to first understand how changes 
in life history parameters will affect the equilibrium density of the different 
classes in the model system (Figured]). Increases in the harem size, h, or 
the overall birth rate leads to an increase the total density of males and 
females (equally since the sex ratio is 50:50). If the male death rate is re- 
duced (relative to the female death rate) then there is an increase in the 
overall density of males (through an increase in infected males) and a de- 
crease in female density. As the male death rate is increased then there is an 
increase in the overall density of females (through an increase in susceptible 
females) and a decrease in male density. Note also that the prevalence of 
infection decreases as the male death rate increases (relative to the female 
death rate). When there is a reduction in the competition parameter for 
males there is an increase in male density through increases in susceptible 
and infected males and female density in both susceptible and infected classes 
is reduced. When there is an increase in the competition parameter for males 
then overall male density and female infected density decreases while female 
susceptible density increases. The overall susceptible density remains con- 
stant as the competition parameter is varied but the proportion of susceptible 
males/susceptible females decreases as the competition parameter for males 
is increased. Also, the prevalence of infection remains relatively constant 
when the male competition parameter is less than the female parameter but 
the prevalence decreases when the male competition parameter is greater 
than the female parameter. 

3.1. Evolving male characteristics. 

The results when the male characteristics are allowed to evolve against 
fixed female parameters are shown in figure [21 As harem size decreases the 
level of resistance to disease in males decreases (the singular value of /3m 
increases). Decreased disease resistance in males also evolves as the male 
death rate increases and the male competition parameter increases. This 
decrease in disease resistance in males is a response to decreased levels of 



prevalence of infection (see figured]) which reduces the need to avoid infection 
(as individuals are less likely to become infected). Male-biased parasitism is 
therefore evident under polyandrous mating systems and when males have a 
higher death rate or suffer more severe competition than females. 

3.2. Coevolving male and female characteristics 

When both male and female characteristics are allowed to evolve variation 
to the harem size does not produce a bias between male and female infection 
rates (figure |3^). Both males and females will evolve increased levels of 
resistance as harem size increases (figure [3]d) in response to the associated 
increases in prevalence. 

When the male death rate exceeds that of females then male biased par- 
asitism can result from coevolution (figure H^). Here the increased death 
rate for males means they have on average a shorter lifespan and so increases 
the possibility of dying from natural causes before becoming infected. This 
is reflected in the fact that the proportion of female/male susceptibles in- 
creases as the male death rate increases. Since there are more susceptible 
females than males it implies that females are more likely to be infected and 
therefore males can afford to pay the cost of an increase to the infection rate. 
When male biased parasitism occurs it increases as harem size decreases. As 
harem size decreases the prevalence of infection decreases and the evolved 
level of transmission increases (disease resistance decreases) for both males 
and females (figure 3). This accentuates the relative difference between male 
and female transmission in a multiplicative manner. 

A similar response occurs when male competition exceeds that of females 
(figure Hb) and can again be attributed to changes in the proportion of fe- 
male/male susceptibles as the male competition parameter increases. When 
the male competition parameter is reduced below that of the female parame- 
ter the evolved level of transmission remains relatively constant. This occurs 
as the prevalence of infection also remains relatively constant and therefore 
there is no selection for a change in resistance to infection. 

3.3. Generality of results for other trade-off's 

We have undertaken the above analysis when the level of parasitism, f3, 
is traded-off against the competition parameter, q{q = g{f3)), and where we 
have imposed a trade-off such that a decrease in parasitism rate for females 
results in an increase in the parasitism rate for males (/3m = 9{Pf))- We 
find the results are analogous to those presented above. The changes in 
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the mating system do not produce a bias in parasitism between males and 
females. Male biased parasitism occurs when the males have a short lifespan 
(or where appropriate suffer increased competition) in comparison to females. 
Again the type of mating system can only accentuate this bias rather than 
cause it. 



4. Discussion 

Differences in the rate of parasitism betwee n sexes and in particular mal e 



biased parasitism i s often found in nature (see ISkorping and JensenI ( 120041 ) ; 
Moore and Wilson! fl2002h ; [Zu3 fll990l l2009f ) ) . We examined whether parasite 
bias between the sexes could arise as a result of the mating system or through 
differences in the underlying life history characteristics between males and 
females through natural selection due to the epidemiological feedbacks that 
they cause. As such we are examining the evolutionary ecological implica- 
tions of life-history and mating system in isolation from the role that they 
may play in sexual selection. Male-biased parasitism was selected for when 
males have a shorter lifespan than females or when males were subject to 
greater competition for resources than females (provided the overall level of 
competition was not too low). Female-biased parasitism will evolve if these 
conditions were reversed. Changes to the mating system did not produce 
a bias in parasitism but could accentuate an existing bias. In particular as 
harem size decreases an existing male biased parasitism is increased as a 
result of a decrease in overall prevalence. We therefore predict more male 
biased parasitism when males have shorter lifespans than females in monog- 
amous or polyandrous species. 

Selection for biased rates of parasitism requires there to be underlying 
differences in the life history characteristics of males and females. Male- 
biased m ortality rates have been repo rted in vertebrate and invertebrate 
systems ( iPromislowl . Il992t iRolfla . |2002| ) and it has been shown to have a 
positive correlation with male-biased parasitism ( iMoore and Wilson! . !2002! ). 
It has therefore been suggested that male-biased paras i tism m ay drive the 
increased mortality rates in males (JMoore and Wilsonl . !2002l ). Our study 
indicates that male-biased parasitism may evolve as a consequence of male- 
biased mortality that was o f course recognised as a p ossible interpretation 
of the empirical findings in !Moore and Wilsonl ( 12002! ). We show that the 
increased mortality in males leads to differences in the population dynamics 
with a greater proportion of susceptible females than males. This reduces 
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the likelihood of infection for males and so they can afford to select for higher 
rates of parasitism. Fundamentally our argument is that non-disease causes 
of higher mortality in males per se may select for the observed decrease in 
immune investment. This increase in mortality can come from processes 
such as increased risks from larger range size or fighting over females. By 
examining how natural selection operates through the ecological feedbacks 
we show in general terms that if males are shorter lived they will invest less 
in resistance. 

The mating system, determined by the choice of the harem size, does 
not directly select for differences in parasitism between males and females in 
our models. However, differences in life-history characteristics that select for 
parasitism bias can be accentuated by the mating system. This is because as 
the harem size reduces, total prevalence levels also reduce leading to selec- 
tion for reduced levels of host resistance. The mating system, which in this 
study acts via the harmonic birth function, effects population density and 
the overall level of disease resistance that evolves. If there is less selection 
for resistance then given an existing bias, it becomes accentuated. Taken as 
a whole, our results on the importance of niating systems are very different 
to those expected from sexual s election (IZukl . ll99Cll . l2009l : IStoehr and Kokkol . 
20061 : iMoore and Wilsonl . |2002| ). Generally polygamous species are expected 
to have stronger selection and therefore more male biased parasitism. Our 
results show that for a 50-50 sex ratio the ecological feedbacks that operate 
due to natural selection affect males and females equally. Therefore, the mat- 
ing system per se does not lead to a sexual bias. Furthermore, the effect of 
mating system on accentuating existing biases runs counter to the prediction 
of the sexual selection idea. Monogamous or polyandrous mating systems are 
more likely to accentuate the bias and therefore show male biased parasitism. 

The mating system can have important consequences for the population 
dynamics that are e xhibited and can lead to cornphcated (cycles, chao s) dy- 



namical behaviour ( Lindstrom and Kokkol . Il998l : iMiller et al.l . l2007bl ). Re 



cently, the population dynamical effects of male-biased parasitism for differ- 
ent case mortalities and both monogamous and polygamous mating systems 
( Miller et al.l . l2007bl ) have been examined. The population dynamics exhib- 
ited (point stability, cycles, chaos) did not show clear trends with increasing 
male-biased parasitism and the outcome depended on a complex interac- 
tion between the ho sts mating system, demography and parasite virulence 



(JMiller et al.l . l2007bl ). Here we focus on the situation where the underlying 



population dynamics are at a stable point equilibrium to allow analysis of the 



10 



fitness expression. We also choose the trade-off to ensure the singular point is 
an evolutionary stable attractor. Studies which examine the evolutionary be- 
haviour for non-equilibrium underlying dynamics are rare but report that the 
evolutionary behaviour would not change for trade-offs where an attractor 



is pr edicted under equilibrium conditions ( iWhite et al.l . l2006l : iHoyle et al 



2010l ). We would therefore expect our finding to extend to non-equilibrium 
underlying dynamics. 

Our result s confirm pr e vious general work on the evolution of resistance 



to parasites ( iMiller et al.l . l2007al : iBoots et al.l . |2009| ) in that we find that 



the evolved level of host resistance increases as the average lifespan of the 
host increases. This result is linked to increased levels of prevalence which 
occur as lifespan increases. In our study the prevalence levels also increase 
as the harem size increases or as the level of competition for resources is 
reduced (until the competition for resources is low) (figure H]). As prevalence 
of infection increases it is beneficial to evolve higher levels of resistance in 
an attempt to avoid infection. Throughout we have assumed that there is 
no long-lived immunity after recovery from infection. In principle this may 
have important consequences with circumstances under which longer-lived 
individuals do not invest more in immunity. As such male biased parasitism 
may be less likely in disease with long lasting immunity, but a full theoretical 
analysis has not as yet been carried out. 

This theoretical study examines the evolution of male-biased parasitism 
in the context of the complex epidemiological feedbacks in disease systems. 
A recent paper has also shown the importan ce of epidemiolog i cal fe edbacks 
to the evolution of male biased parasitism ( iResif and Amosl . l2010l ). In a 
comprehensive study, they examine how differential exposure between males 
and females affects various aspects of investment in immun ity under a range 
of tra de-offs including one between recovery and lifespan ( IResif and Amosl . 
20101). The model includes diploid genetics mapped onto a quantitative trait 
and fundamentally includes the epidemiological feed backs caused by differ- 
ent investments under monogamous mating systems. iResif and Amosl (12010! ) 
show that reduced investment in males can evolve when there is more ex- 
posure to parasites (achieved by imposing differences in some of the disease 
characteristics between males and females). Their results further emphasise 
the importance of epidemiological feedbacks. Our study does not impose un- 
derlying differences in disease characteristics but focuses on the role of host 
life history and mating system. We have shown that male-biased parasitism 
can evolve as a consequence of sexual differences in life history characteristics 
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that produce a greater proportion of susceptible females than males. Our re- 
sults extend to different choices of trade-offs. Future studies should extend 
the analysis to examine the importance of the choice of underlying infectious 
disease framework and the representation of the two-sex birth function that 
may include assessing the effects of a non-equal sex ratio. Throughout we 
are focussing on the role of natural selection in the context of epidemiological 
feedbacks. Future work could combine this approach with models of sexual 
selection in order to gain a full understanding of the mechanisms that un- 
derpin m ale biased parasi t ism. T he combined genetic and quantitative trait 
model of iResif and Amod (I2OIOI ) could be extended to provide a framework 



in which to examine these different processes. 
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Figure 1: Changes in the density of the different classes as other parameters are varied 
where "o" represents S*/, "•" represents //, "n" represents Sm and "■" represents Im- 
Parameter values (unless varied in the figure) are Cm ~ Cf — 3.79, Qm = Qf = 0.25, 
bm ~ bf — 1, am — oif — 1, j3m — Pf — 1-6625, h — 1. In a) ft, a is varied, b) c™ is varied, 
c) bm is varied and d) qm is varied. In a) and b) Sf — Sm and If = Im so the results are 
just shown for males. 
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Figure 2: Change in the singular value of disease transmission for males, /3^, when males 
can evolve ( and females do not evolve). Parameters values (unless varied in the figure) 
are c/ = 3.79, qm, — Qf — 0.25, bm ^ bj = 1, am — Oif — ^t Pf = 1.6625, h = 1. The 
parameters for the trade-off (equation [3]) are [Smin = 1, l^max — 20, Cmin — 2, Cmax = 5 
and 7 = 40. In a) h is varied, b) 6™ is varied and c) q™ is varied. 
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Figure 3: Change in the coevolutionary singular value of disease transmission, /3*„ and /3i 
plotted against changes in harem size. Parameter values are Qm — q/ — 0-25, &,„ — bf — 1, 
ttm = «/ = !. The parameters for the trade-off (equation [S]) are as in Figl2] a) shows the 
relative values of (3m and /?/ at the coevolutionary singular point and b) shows the actual 
value of /3*j at the coevolutionary singular point (note the value of /3i is identical here). 
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Figure 4: Change in the coevolutionary singular value of disease transmission /3*j and 
/3i plotted against changes in a) male birth rate and b) male susceptibility to crowding. 
Parameter values (when not varied in the figure) are bm = bf = 1, Qm = Qf = 0.25, 
ctm = ctf = ^, Pm = /3/ = 1.6625. The parameters for the trade-off (equation [3]) are as in 
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